Signatures of the Tricritical Point in QCD 
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Several approaches to QCD with two massless quarks at finite temperature T and baryon chemical 
potential /i suggest the existence of a tricritical point on the boundary of the phase with sponta- 
neously broken chiral symmetry. In QCD with massive quarks there is then a critical point at the 
end of a first order transition line. We discuss possible experimental signatures of this point, which 
provide information about its location and properties. We propose a combination of event-by-event 
observables, including suppressed fluctuations in T and /i and, simultaneously, enhanced fluctuations 
in the multiplicity of soft pions. 



In QCD with two massless quarks, a spontaneously 
broken chiral symmetry is restored at finite temperature. 
It can be argued that this phase transition is likely 
second order and belongs to the universality class of 0(4) 
spin models in 3 dimensions. If this transition is indeed 
second order, QCD with two quarks of nonzero mass has 
only a smooth crossover as a function of T. Although 
not yet firmly established, this picture is consistent with 
present lattice simulations and many models. 

At zero T several models suggest that the chiral 
symmetry restoration transition at finite /i is first order. 
Assuming that this is the case in QCD, one can easily 
argue that there is a tricritical point in the T/z phase 
diagram, where the transition changes from first to sec- 
ond order. The nature of this point can be understood 
by considering the Landau- Ginzburg effective potential 
for the order parameter of chiral symmetry breaking, 
(j) = (cr,7r) - {Tpil)): 
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The coefficients a, 5, and c > arc functions of /i and 
T . The second order phase transition line described by 
a = at 6 > becomes first order when b changes sign. 
The critical properties of this point can be inferred from 
universality |3,0, and the exponents are as in the mean 
field theory (1). 

In real QCD with nonzero quark masses the second 
order phase transition becomes a crossover and the tri- 
critical point becomes a critical (second order) end-point 
of a first order phase transition line. Universality argu- 
ments [^,0 also predict that the end-point E in QCD with 
small quark masses is shifted with respect to the tricrit- 
ical point P towards larger /i as shown in Fig. ^ It can 
also be argued |^,^ that the point E is in the universality 
class of the Ising model in 3 dimensions, because the a 
is the only field which becomes massless at this point. 
(The pions remain massive because of the explicit chiral 
symmetry breaking by quark masses.) In this paper we 
discuss experimental signatures of this critical end-point. 

The position of the points P and E in two-flavor QCD 
was estimated recently using two different models (a 
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FIG. 1. The schematic phase diagram of QCD. The dashed 
lines represent the boundary of the phase with spontaneously 
broken chiral symmetry in QCD with 2 massless quarks. The 
point P is tricritical. The solid line with critical end-point E 
is the line of first order transitions in QCD with 2 quarks of 
small mass. The point M is the end-point of the nuclear liq- 
uid-gas transition probed in multifragmentation experiments. 
The superconducting phase of QCD Pj^,^, marked SC, is not 
relevant to our discussion. 



Nambu-Jona-Lasinio model respecting the global sym- 
metries of QCD and a random matrix model 0) as 
Tp - 100 MeV and /xp ~ 600 - 700 MeV. These are only 
crude estimates, since they are based on modeling the 
dynamics of chiral symmetry breaking only. 

The third (strange) quark has an important effect on 
the position of the point P and, therefore, of the point E. 
At /i = 0, if the strange quark mass rus is less than some 
critical value rriss, the second order finite T transition 
becomes first order. This leads to a tricritical point in 
the Trus plane @,|l]j|. Theoretically, the origin of this 
point is similar to the one we are discussing. In terms of 
eq. (|l]) the effect of decreasing is similar to the effect 
of increasing /i: the coefficient b becomes negative. What 
is important is that, unlike m^, /i is a parameter which 
can be experimentally varied. 

Clearly, the physics of the T/i plane is as in Fig. |l| 
only for mj > 77153. For 77is < 77753, the transition is 
first order already at /x = 0, and, presumably, remains 
first order at all nonzero /x ||l^ . As 7775 is reduced from 
infinity, the tricritical point P of Fig. |l| moves to lower 
/i until, at 777s = '77s3, it reaches the T-axis and can be 
identified with the tricritical point in the r777s plane. The 



1 



two tricritical points are continuously connected. We as- 
sume that iTis > rriss which is consistent with the lattice 
studies of ref. . What is important for us is that the 
qualitative effect of the strange quark is to reduce the 
value oi Hp, and thus of fiE, compared to that in two- 
flavor QCD, since /ip = at nis = rriss- This shift may 
be significant, since lattice studies show that the physical 
value of rris is of the order of TOss. 

Analysis of particle abundance ratios in central heavy 
ion collisions indicates that chemical freeze-out hap- 
pens near the phase boundary, at a chemical potential 

- 500 - 600 MeV at the AGS (11 GeV-A), while at 
the SPS (160 - 200 GeV-A) it occurs at a significantly 
lower /X ^ 200 MeV. In view of the effect of the strange 
quark just discussed, the estimated position of P and E 
should be shifted from ~ 600-700 MeV to lower 
11. Thus, it may well be between the SPS and the AGS 
values of ^, and therefore the point E may be accessible 
at lower energy or non-central collisions at the SPS. 

The strategy for finding the point E which we propose 
is based on the fact that this is a genuine critical point. 
Such a point is characterized by enhanced long wave- 
length fluctuations which lead to singularities in all ther- 
modynamic observables. In the liquid-gas phase transi- 
tion in water, critical opalescence signals the universal 
physics unique to the vicinity of the critical point. The 
signatures we propose can play an analogous role in QCD. 

It is important to have control parameters which can 
be adjusted to vary the /i at which the system crosses 
the transition region, as shown in Fig. 2. For exam- 
ple, increasing the energy of the collision decreases this 
fi. A somewhat similar effect may be achieved by in- 
creasing centrality. A third possibility would be to slice 
each event in rapidity, since /i will be greater at higher 
rapidity. This strategy could be useful at RHIC, if E 
were to lie at lower ^ than is accessible at the SPS. We 
will call the control parameter which is varied "x" , and 
take increasing x to mean increasing collision energy or 
centrality or decreasing rapidity. Scanning in centrality 
will almost certainly be the easiest, since in any given 
run events with all impact parameters are present. How- 
ever, scanning in energy yields a large variation in the fi 
at which the transition is crossed, whereas scanning in 
centrality may only provide fine tuning. 

In this work we do not discuss initial equilibration and 
we choose to define the initial point, l{x), as the point at 
which compression has ended, most of the entropy is al- 
ready produced, and approximately adiabatic expansion 
begins. The system will then follow some trajectory in 
the Tfi plane characterized by the ratio of the baryon 
charge density to the entropy density, n/s, which is (ap- 
proximately) conserved. Three trajectories arc shown 
schematically in Fig. ||. (For realistic hydrodynamical 
calculations and discussion see, e.g., refs. |15 
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Recall that the first order line in the T/i plane is actu- 
ally a whole region of mixed phase, with the hidden pa- 
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FIG. 2. Schematic examples of three possible trajectories 
for three values of x on the phase diagram of QCD (see. Fig. 
|l]). The points I, S, H and F on different trajectories are 
marked with different symbols. The dashed lines show the 
locations of the initial, I, and final, F, points as x is increased 
in the direction shown by the arrows. 



rameter being the volume fraction of the two coexisting 
phases. The zig-zag shape occurs because the trajectory 
exits the mixed phase region at a point with the same 
value of the conserved n/s that it had upon entering. 
Because n/s is discontinuous at the first order line, this 
requires increasing T and decreasing /i as latent heat is 
released p5| , p^ . In Fig. ^, we use the following notation: 
S{x) for the "softest" point, II(a;) for the "hottest" point 
and F(a;) for the final thermal freeze-out after which no 
scattering occurs. (Note that at small values of x, at 
which the transition is first order, the trajectories are, 
in fact, likely to begin within the mixed phase region. 
The special case when 1(2;) coincides with S(a;) leads to a 
local maximum of the QGP hfetimc HIH- which may 
be important for J/ip suppression [p^.) Increasing x will 
yield trajectories shifted to the left in Fig. ||, traversing 
the transition region at lower fi and higher T. 

The existence of the end-point singularity, E, leads to 
the phenomenon which we refer to as the "focusing" of 
trajectories towards E. The initial point l{x) and the be- 
ginning of the zig-zag S(a::) depend on the control parame- 
ter X more strongly than the zig-zag end-point II(a;). The 
reason for this is that the point H(a;) is always closer to 
E than S{x) (see Fig. |^). This focusing effect implies 
that exploring physics in the vicinity of the end-point 
singularity may not require a fine-tuned x. This situa- 
tion resembles that in low energy nuclear collisions, in 
which the first order liquid-gas phase transition also has 
a critical end-point at a temperature of order 10 MeV 
po| , pT| (point M on Figs. |^ and|2|). In such experiments, 
one varies control parameters to maximize the probabil- 
ity of multi-fragmentation. It was noticed long ago l2C| ] 
how surprisingly easy it is to hit the critical region from 
a wide range of initial parameters. 

Another aspect of the "focusing" arises via the diver- 
gence of susceptibilities, such as the specific heat capac- 
ity cv = Tds/dT, at the endpoint E. As a result, the 
trajectories which pass near the critical point will linger 
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there longer. This makes it likely that final freeze-out 
occurs at a temperature quite close to Tg, rather than 
below it. So, while scanning in some control parame- 
ter X and measuring the positions of the points F(a;), we 
may expect to find a bump in the vicinity of the point 
E. (See the lower dashed curve on Fig. ^) At this point 
it is instructive to consider the dependence on another 
control parameter, the atomic weight A of the colliding 
nuclei. If A were infinite, the point F would be close to 
Tp = 0, /ii? = mAT. Thus, for A large enough, the dot- 
ted curve in Fig. ^ moves down and the bump, and also 
all the other signatures described below, fade away. Ex- 
perimentally, the A dependence of the point F has been 
established recently by the analysis of flow [|l6| , Coulomb 
effects and pion interferometry For example, in 
central S+S coUisions at SPS Tp « 140-150 MeV, while 
for Pb-KPb it is only Tp k. 120MeV. 

We shall now discuss the signatures which directly re- 
flect thermodynamic properties of the system near its 
critical point and are not very sensitive to the details of 
the evolution. With the advent of wide- solid- angle de- 
tectors like NA49 at CERN, it is now possible to make 
event-by-event measurements of observables which are 
proxies for the freeze-out T and /i [Q. We argue that 
the event-by-event fluctuations in both quantities should 
be anomalously small for values of x such that the sys- 
tem passes near the critical point. As has been suggested 
p5[ , event-by-event fluctuations of T can be related by 
basic thermodynamics to the heat capacity at freeze-out 

The quantity Cy is extensive, so AT ^ 1/VN as ex- 
pected, where N is the number of particles in the system. 
If the speciflc heat cy diverges, the coefficient of 1/Vn 
vanishes and fluctuations of T are suppressed. For freeze- 
out in the crossover region, or in the hadronic phase just 
below the first order transition, cy is finite. (If freeze-out 
were to occur from the mixed phase, some linear combi- 
nation of the two susceptibilities would be relevant.) As 
the critical point is approached from either the left or the 
right, Cy diverges and AT^/N decreases. Other suscep- 
tibilities, in particular, —d'^il/dfi^, are also divergent. 
This implies that fiuctuations of /i are also suppressed 
at the critical point. Experimentally, AT can be found 
via event-by-event analysis of pT spectra [p5|| . Fluctu- 
ations in fi correspond to event-by-event fluctuations in 
the baryon-number-to-pion ratio. The fluctuations in any 
experimental observable will receive contributions in ad- 
dition to the thermodynamic ones we describe, for exam- 
ple from fluctuations in the flow velocity. We therefore 
expect that as the collision energy is increased so that the 
freeze-out point moves from right to left past the critical 
point, we will flnd minima (but not zeroes) in the widths 
of the distributions of those event-by-event observables 
which are well-correlated with T and /k. 



Using universality, we can predict the exponents for the 
divergent susceptibilities at the point E. Very naively, one 
might think that the exponent describing the divergence 
of Cy is a, which is small for the 3-dimensional Ising 
model universality class: a ^ 0.12. In fact the exponent 
for Cy is significantly larger. This and the exponent for 
the /i-susceptibility are determined by finding two direc- 
tions, temperature-like and magnetic-field-like, in the Tfi 
plane near point E, following the standard procedure for 
mapping a liquid-gas transition onto the Ising model [ p6[ . 
The two linear combinations of T — Tp and ji — iip corre- 
sponding to these directions should then be identified (in 
the sense of the universality) with the temperature, or 
t = T — Tc, and the magnetic field, h, in the Ising model. 
One can easily understand this by realizing that the t- 
like direction should be tangential to the first-order line 
at the point E. Then Cy and —d^il/dfi^ are different lin- 
ear combinations of the t-\ike and ft,-like susceptibilities. 
In both linear combinations, the divergence of the /i-like 
susceptibility will dominate because 7 « 1.2 ^ a « 0.12. 
The exponent for the divergence of the h-\ike suscepti- 
bility as a function of the distance, £, from the point E 
will depend on the direction along which one approaches 
this point. For almost all directions it will be given by 
7//3(5 « 0.8 (except for exactly the i-like direction, where 
it is 7). As a result, for points on the T/i plane along a 
generic line through E one finds 

{ATf ^ (Ap.)^ ^ £"-^ (3) 

sufficiently close to E. Therefore, the fluctuations of T 
and n are considerably suppressed when the freeze-out 
occurs near the critical point. 

We turn now to direct signatures of the long- 
wavelength fluctuations of the massless a field. For the 
choices of control parameters x such that freeze-out oc- 
curs at (or near) the point E, the cr-meson is the most nu- 
merous species at freeze-out, because it is (nearly) mass- 
less and so the equilibrium occupation number of the 
long-wavelength modes [T/lo) is large. Because the pi- 
ous are massive at the critical point E, the cr's cannot 
immediately decay into tttt. Instead, they persist as the 
density of the system further decreases. It is important 
to realize that after freeze-out, one can (by definition) 
approximately neglect collisions between particles. Col- 
lective effects related to forward scattering amplitudes 
cannot be neglected. That is, although the particles no 
longer scatter, their dispersion relations will not be given 
by those in vacuum until the density is further reduced 
by continued expansion. 

During the expansion, the in-medium sigma mass rises 
towards its vacuum value and eventually exceeds the tttt 
threshold. As the cttttt coupling is large, the decay pro- 
ceeds rapidly. This yields a population of pious with 
small transverse momentum, pT < TO,r- Because this 
process occurs after freeze-out, the pions generated by it 



3 



do not get a chance to thermalize. Thus, the resulting 
pion spectrum should have a non-thermal enhancement 
at low pt which is largest for freeze-out near E where the 
(t's are most numerous. 

We now propose a specific signature of the endpoint 
visible in the fluctuations of the pions resulting from 
the (formerly) long wavelength modes of the a field. 
For freeze-out close enough to E that the sigma mass 
at freeze-out is less than T, the thermal fluctuations of 
the number, N„, of a particles are determined by the 
classical statistics of the field cr, rather than by Poisson 
statistics of particles. Therefore, (iV^) - {N„Y ^ i^a)^, 
rather than (N^). Thus, we expect large event-by-event 
fluctuations in the multiplicity and distributions of the 
soft pions: N^^ w 2iVo-. Due to critical slowing down, 
non-equilibrium effects may further enhance these fluc- 
tuations. Thus, these pions could be detected either di- 
rectly as an excess in the p^-spectra at low px, or via 
increased event-by-event fluctuations at low px, or by 
an increase in Hanbury-Brown-Twiss (HBT) correlations 
due to the larger number of pions per phase space cell at 
lowpT- 

To conclude, we propose that by varying control pa- 
rameters such as the collision energy and centrality, one 
may find a window of parameters for which the Tfi tra- 
jectories pass close to the critical point E. Enhanced crit- 
ical fluctuations of the a field and the associated thermo- 
dynamic singularities lead directly to the signatures we 
propose. When the freeze-out occurs near the point E, 
we predict large non-thermal multiplicity and enhanced 
event-by-event fluctuations of the soft pions. In contrast, 
the event-by-event fluctuations in both T and /i, as de- 
termined using pions with px JSJ TOtt , will be anomalously 
suppressed. Both effects should disappear if the atomic 
weight A is very large. No one of these signatures is 
distinctive in isolation and without varying control pa- 
rameters. Several of them seen together and seen to turn 
on and then turn off again as a control parameter is var- 
ied monotonously would constitute a decisive detection 
of the critical point. 

What would we learn about QCD if such a point is 
found? First, we would learn that there is a genuine crit- 
ical point in the Tfj, plane in nature. Second, we would 
learn that m^. > m^a in nature, and the /i = thermal 
transition is a crossover for physical quark masses, rather 
than a first-order phase transition. Third, the experimen- 
tal discovery of the critical end-point E would mean that 
if the light quark masses were set to zero, there would be 
a tricritical point P in the phase diagram of QCD. 
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